Two-Dimensional Electron Gas in InGaAs/InAlAs Quantum Wells 
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We designed and performed low temperature DC transport characterization studies on two- 
dimensional electron gases confined in lattice-matched Ino.53Gao.47As/Ino.52 Alo.4sAs quantum wells 
grown by molecular beam epitaxy on InP substrates. The nearly constant mobility for samples 
with the setback distance larger than 50nm and the similarity between the quantum and trans- 
port life-time suggest that the main scattering mechanism is due to short range scattering, such 
as alloy scattering, with a scattering rate of 2.2 ps . We also obtain the Fermi level at the 
Ino.53Gao.47As/Ino.52 Alo. 48 As surface to be 0.36eV above the conduction band, when fitting our 
experimental densities with a Poisson-Schrodinger model. 



Two dimensional electron gas (2DEG) confined in 
Ino.53Gao.47As in lattice matched InGaAs/InAlAs/InP 
heterostructures and superlattices appear in many tech- 
nologically important areas ranging from high speed 
electronics Q, optoelectronics 0, an d spintronics0, 
It is also an attractive 2DEG system for the study 
of disorder induced quantum phase transitions @, 0, 
H, U E3, [n|. While there have been several earlier 
works characterizing electronic properties of 2DEGs in 
Inn . 53 Gao . 47 As/Inn . 52 Alo .48 As hetero j unctions [12], [TJ, [bj, 
Ho| . Inp.sa Gao 47A . S/I1 1P heterojunctions and quantum 
wells (QW) [la. ll7L Hq | . there were few systematic studies 
characterizing 2DEGs in Ino.53Gao.47As/InpR2Alo.48As 
QWs. Since many modern structures 0,0 HIM are now 
based on Ino.53Gao.47As/Ino.52Alo.4sAs QWs, such char- 
acterization is of fundamental and technological interest. 

In this letter we report the characterization of elec- 
tronic properties of 2DEG in a series of lattice matched 
InGaAs/InAlAs QWs grown by molecular beam epitaxy 
(MBE) on InP substrate (here and after in our paper we 
use abbreviations InGaAs for Ino.53Gao.47As and InAlAs 
for Ino. 52 Alo. 48 As). Systematic investigations of 12 such 
wafers with varying design parameters in the doping lay- 
ers have yielded important information not only about 
carrier mobility and scattering, but also about how dop- 
ing determines the carrier densities, from which we were 
also able to determine the location of the Fermi level at 
the InGaAs surface. 

The schematics of the samples is depicted in Fig. ^a) 



and the parameters for each sample are summarized in 
Table H The 2DEG resides in a 20nm-wide InGaAs QW. 
Two Si (S-doped layers are placed in the InAlAs barrier 
to one side (closer to the surface) of the QW. The three 
design parameters that were varied are the doping den- 
sities (N t and Nb) in the top and bottom dopant layers 
respectively, and the distance d from the bottom dopant 
layer to the (top) edge of the QW. We fabricated stan- 
dard Hall bars with Indium ohmic contacts. We tried to 
measure all the samples at dark. Except a few samples 
(6 — 9) most of them need to be illuminated to create 
a 2DEG. For these samples we illuminated for sufficient 
time with an LED to create a 2DEG with the highest 
possible mobility. We measured the magnetoresistance 
R xx and the Hall resistance R xy as a function of the 
perpendicular magnetic field (B) for different tempera- 
tures. The results are shown in Fig. ^ for the sample 5 
(see Table Similar results were obtained for the other 
samples. From the measured data we obtain the areal 
density n2D and the mobility /1 of the electrons, from 
which we can extract the transport lifetime r t — fim* /e. 
We used an effective mass m* for Ino.53Gao.47As of 0.043 
times the bare electron mass[l9j. From the onset of 
SdH oscillations we extracted also the quantum lifetime 
(r q ) by using a Dingle style analysis. The amplitude 
(AR) of the envelope function of the SdH oscillations 
was found to be well described by the conventional Ando 
formula [HEl Ai?sinh(A T )/4 J R ^T = e-*/ UcT «, where 
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FIG. 1: The main panel shows p xx and pnaii as a function 
of the magnetic field for temperatures from 22 (red trace) 
to 1300 mK (yellow trace), for the sample 5. The inset (a) 
shows the schematic diagram of the samples. The substrate 
is semi-insulating InP. The inset (b) illustrate how we extract 
the amplitude of the SdH oscillations as a function of B, at 
1300 mK. Finally we plot the amplitude in a Dingle plot (c) 
to obtain the scattering quantum time r q . 



TABLE I: Sample parameters, electron densities U2D, mobil- 
ities fj,, and the transport r t and quantum r q scattering times 
for 12 different structures at 4.2 K. 



Sample 


d 


N t 


N b 


ri2D 




Tt T q 

(ps) (ps) 


number 


(nm) 


(lO 11 /™ 2 ) 


(10 11 /cm 2 ) 




(em 2 /V») 


1 





2 


1 


3.1 


4500 


0.11 0.17 


2 


20 


2 


1 


3.0 


13000 


0.31 0.22 


3 


50 


2 


1 


1.9 


16000 


0.39 0.31 


4 


50 


1 


1 


1.7 


14500 


0.35 0.31 


5 


50 


2 


1 


1.9 


16000 


0.39 0.32 


6 


50 


5 


1 


2.2 


15500 


0.38 0.31 


7 


50 


10 


1 


3.3 


15000 


0.37 0.31 


8 


50 


10 





2.1 


15000 


0.37 0.31 


9 


50 


5 


0.5 


1.7 


15000 


0.37 0.31 


10 


150 


2 


1 


1.6 


15500 


0.38 0.33 


11 


200 


10 


1 


1.7 


15000 


0.36 0.31 


12 


300 


10 


1 


1.5 


15500 


0.38 0.31 



quency and i?o represents the resistance at zero applied 
magnetic field, for a given value of temperature. To ob- 
tain the amplitudes ( AR) , for each temperature we fitted 
the envelope of the SdH oscillations to a pair of polyno- 
mials as showed in Fig. ^b) for T=1300 mK. Next we 
obtained AR just subtracting both polynomials and we 
plotted Ai?-sinh(AT)/4i?o^T versus 1/B in alog-x graph 
as in Fig.^c). From the previous mentioned Ando for- 
mula, we performed a linear fit to achieve the slope (s) of 
the Dingle plot [Fig. ^c)] and obtain the quantum life- 
time for the given temperature as r q — — (it ■ m*)/(e • s). 
Table [I] summarizes our results measured at 4.2 K, to- 
gether with relevant parameters of the samples. 

Our measurements show no indication of parallel con- 
duction nor any presence of a second subband. We 
have also confirmed this by solving self-consistently the 
Schrodinger and Poisson equations to calculate the sub- 
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FIG. 2: The density of the 2DEG (nan) is plotted for the dif- 
ferent samples. Lines are guides to the eye. We used Vs=0.36 
(eV) as the best fit to the experimental data for the different 
samples. 



band energy levels. The calculation shows that only one 
subband is occupied and the second subband is more than 
40 meV above the Fermi level. To understand the origin 
of the electrons forming the 2DEG, we solved analyti- 
cally the electrostatic Poisson equation for our structure 
following similar procedures as in Ref. |^] and assuming 
full ionization of dopants. This leads to the following 
analytical expression for n 2 o for our structure, 



n 2D 



N b {c + h+ l 2 ) + N t (c + l 2 ) + {V B e e B ) 

C + h+ln+d + Z ' 



(1) 



where d, N t , N b are as defined earlier, c=50nm is the 
thickness of the InGaAs cap layer, i2=18nm is the dis- 
tance from the top Si <5-doped layer to the upper-end 
of the InAlAs and Zi=150nm the distance between the 
two Si (5-doped layers; 6b=14.2 is the relative dielectric 
constant of InGaAs 22] and V B is the offset of the fermi 
level at the InGaAs surface with respect to the conduc- 
tion band edge. In the denominator, z is an effective 
"quantum" depth of the 2DEG[2lJ which turns out to be 
approximately the QW width (20nm) in our case (the de- 
viation is negligible compared to the contribution from 
other terms in the denominator). If we plug- in all the 
relevant numerical values, we get the following dimen- 
sionless formula for n 2 o'- 



n 2D 



173 N b + 23N t + 785V B 
d+ 193 ' 



(2) 



where n 2 D, Nf,, and Nt are in units of 10 11 cm -2 , d in nm, 
and Vb in eV. The only free parameter in our model, V B , 
is the Fermi level relative to the conduction band edge 
(CBE) at the InGaAs surface. Due to surface states, 
the local Fermi level at a semiconductor surface is often 
pinned regardless of doping and carrier density. For ex- 
ample, the surface Fermi level for GaAs is about ~ 0.7 eV 
below the CBE. For InGaAs, such information is largely 
unknown, hence we have fitted this parameter using our 
n 2 D data and the best fit is obtained for V B — 0.36 eV. 
Here, the surface Fermi level is above the CBE, which 
is very similar to the InAs case [23|. except that in our 
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samples, it appears that the surface carriers are not mo- 
bile enough and do not contribute significantly to the 
transport. 

In Fig. |2 n2D is plotted as a function of sample num- 
ber, with Vb = 0.36eV. Clearly, there is an excellent 
agreement, between this fit (J5J and the experimental 
data. This also confirms that n^D is only a function 
of d, Nb and N t and indicates that unintentional dop- 
ing from the residual or background impurities does not 
appear to be significant in our samples. Indeed, intro- 
ducing background impurities in our model would lead 
to a decrease of the features seen in Fig. [21 and therefore 
does not fit our data as well. From the data in Table [I] 
we can notice that for d larger than 50nm, the mobilities 
and the lifetimes arc independent of the doping param- 
eters (d, Nb and Nt). In Fig. [3| we plot the quantum 
scattering rate r" 1 (triangles) and the transport scatter- 
ing rate T t -1 (circles), measured at 4.2K for samples 1-3 
and 10-12, as a function of d, the distance from the bot- 
tom doping layer to the quantum well. Both scattering 
rates (r _1 ) shows a fast decrease at small values of d 
(below ~30nm), indicating that the dopants provide ef- 
ficient scattering for electrons at short d. For d of 50 nm 
or more, however, t _1 becomes independent of d. From 
samples with G?=50nm but different Nt and Nb, we found 
t _1 to be also independent of the doping densities, as 
is shown, for example, in the inset (a) of Fig. [3] where 
t _1 is plotted against Nt at d=50 nm. Our findings 
indicate that, for d >50nm, the dopant layers are not 
the major source of carrier scattering in these structures. 
The common and reproducible value of mobility (^15000 
cm 2 /Vs) observed for wafers (with d >50 nm) from dif- 
ferent MBE growth suggests that the mobility, or carrier 
scattering originates from some intrinsic, non-doping re- 
lated scattering in our structures. Moreover, since the 
dependence of the mobility on the 2DEG density is very 
weak as seen in Table we believe that an important 
source of scattering is due to the random alloy scatter- 
ing potential^ which is expected to have a weak density 
dependence 24]. Such intrinsic alloy scattering is clearly 
very important in the InQ.53Gao.47As channel, as claimed 
previously to dominate the low temperature scattering in 
InGaAs/InAlAs heteroiunctions jla yjj . We further be- 
lieve that surface roughness is less important since these 
structures are lattice matched MBE grown and because 
surface roughness would lead to a stronger dependence of 
mobility on density. Hence the main source of disorder 
is short-ranged in contrast to charged doping disorder, 
which is long ranged, in relation to the Fermi wavelength. 
The short range nature of the dominant scattering mech- 
anism for our samples with d >50nm is consistent with 
our observation that the quantum life time is similar to 
the transport lifetime. We have measured the scattering 
rate dependence at lower temperatures (T), as shown in 
inset (b) of Fig. [3| for representative data in sample 5, 
from which we extract a low T limiting value of the scat- 
tering rate to be ~2.2 ps _1 . 

Since most samples have to be illuminated in or- 
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FIG. 3: The dependence of the quantum (triangles) and trans- 
port (circles) scattering rates (t _1 ) with the distance d of the 
doping layer to the quantum well . The insets shows the de- 
pendence of t -1 with (a) the amount of doping Nt and (b) 
temperature for the case of d=50 nm. 



der to obtain the optimum mobility, it is likely that 
DX centers play a role |25j, |26|, |27j. We therefore es- 
timated their contribution by fitting, for the sample 
5, the temperature dependence of the total density to 
n 2 D = n free + n DX where rif ree = 1.17 x 10 cm and 
n DX {T) = n' DX exp [{E DX -E F )/kT]. The best fit is 
obtained for n' DX = 2.5xl0 10 cm -2 andE DX = Ef— 21.5 
meV (below the Fermi energy). Hence, DX centers could 
explain the observed increasing of both the 2DEG den- 
sity and of the scattering rate with increasing tempera- 
ture. Indeed, at higher temperatures, more carriers are 
activated, which will also leave the DX centers unsatu- 
rated and lead to the increased scattering rate as seen in 
Fig. lib). 

In conclusion, we have studied the two- 
dimensional electron system confined in MBE-grown 
Ino.53Gao.47As/Ino.52Alo.4sAs quantum wells. We have 
measured 12 different wafers specifically designed and 
grown to investigate transport properties of this type of 
material. We have obtained an analytical expression for 
the ri2D of our samples with a very good agreement with 
the measured values, showing that uid depends only 
on parameters of the Si <5-doped layers (d, N t and Nb). 
We obtain an excellent fit to our experimental densities 
assuming the Ino.53Gao.47As surface Fermi energy to be 
0.36 eV above the conduction band. For setback distance 
d of 50nm or more, quantum and transport scattering 
rates are independent of parameters of the dopant layers 
and are likely mainly due to short-range scattering, such 
as alloy scattering, for which a scattering rate of 2.2 
ps _1 is extracted at 22 mK. 
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